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RESUMEN 

Favor de proporcionar un resumen en espaiiol. If you are unable to translate your abstract into 
Spanish, the editors will do it for you. I review the effects of magnetic fields on star formation in galaxies. 
This includes the effects of the magnetorotational instability (MRI) at galactic scales, magneto- Jeans and 
swing instabilities, Parker instabilities, and the effects of magnetic fields on the evolution of supernova-driven 
turbulence. I argue that currently turbulent support by the MRI appears likely to be the most important of 
these processes to regulating star formation. 

ABSTRACT 

I review the effects of magnetic fields on star formation in galaxies. This includes the effects of the magnetoro- 
tational instability (MRI) at galactic scales, magneto-Jeans and swing instabilities, Parker instabilities, and 
the effects of magnetic fields on the evolution of supernova-driven turbulence. I argue that currently turbulent 
support by the MRI appears likely to be the most important of these processes to regulating star formation. 
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1. GLOBAL STAR FORMATION 

The question of what regulates star formation in 
galaxies remains unsolved, with at least four main 
scenarios, and a number of variations on some of 
them. To understand the role of magnetic fields, we 
must first briefly review these ideas. 

The first is that global star formation is primarily 
controlled by gravitational instability of the avail- 
able gas in the disk. Whether this instability is 
global, controlled by the combined potential of the 
gas and stars Rafikov (2001), or local, controlled 
by the behavior of individual GMCs, is still argued. 
The basic idea of global instability was described by 
Elmegreen (2002), and supported by numerical mod- 
els (Kravtsov 2003; Li et al. 2005a). Observational 
evidence for a direct correlation between star forma- 
tion rate and gravitational instability in the LMC 
is given by Yang et al. (2007). Local instability in 
molecular clouds has been argued to be the rate lim- 
iting step by Krumholz (2005). 

A second idea recently elucidated by Shu et al. 
(2007) is that magnetic regulation dominates star 
formation. The results of Kim & Ostrikcr (2006) and 
Shetty & Ostriker (2006) establish that the magnetic 
field can control the morphology of star formation 
in spiral galaxies. The suggestion is then made by 
Shu et al. (2007) that the rate limiting step in star 
formation is the accumulation of sufficient mass to 
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form a supercritical cloud. 

Star formation controlled by a threshold column 
density has been advocated by Schaye (2004). He 
argues that gas above the threshold will be able to 
cool quickly by forming molecules, which will drive 
gravitational instability. However, this would pre- 
dict uniformly high molecular fractions at the crit- 
ical column density, contrary to the observations of 
Martin & Kennicutt (2001). 

Finally, self-regulated star formation remains on 
the table (e.g. Silk 1997). In this scenario, the su- 
pernovas produced by star forming regions determine 
the level of turbulence in the surrounding gas (e.g., 
Slyz et al. 2005), which determines its level of grav- 
itational instability and thus how many stars it can 
form. Less star formation lowers the level of tur- 
bulence, which in turn increases the star formation 
rate. However, this cannot explain the lack of star 
formation far from star forming regions and massive 
stars, e.g. in the outer disks of normal galaxies, or in 
low surface brightness galaxies. 

2. GRAVITATIONAL INSTABILITY 

Global disk instability models postulate that star 
formation happens wherever gravitational instabil- 
ity (Gammie 1992; Rafikov 2001) of the combined 
collisionless stars (Toomre 1964) and collisional gas 
(Goldreich & Lynden-BcU 1965) in the disk sets in 
during the collapse of disk galaxies. For a gas disk, 
the criterion for instability is 
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where k is the epicychc frequency, Cg the speed of 
sound, and Sg the surface density of the gas disk. 
For a coUisionless stellar disk, the equivalent crite- 
rion is 

^ ^ < 1.07, (2) 

where as is the stellar velocity dispersion in the ra- 
dial direction, and is the surface density of the 
stellar disk. [Note that following Rafikov (2001) 
we use a factor of tt in the definition of Qs rather 
than 3.36, shifting the instability criterion to slightly 
higher value.] Define the dimcnsionlcss quantities 
q = kas/n and R = Cg/cTs- Then the instability cri- 
terion for the combined disk of gas and stars is given 
by 
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where Iq is the Bessel function of order 0. The 
central determining factor in global disk instability 
models is then the magnitude of the turbulent sup- 
port of gas. 

This support cannot come solely from star for- 
mation, however. The 21 cm observations of Petric 
& Rupen (2007) represent the highest-resolution ob- 
servation of velocity dispersion across a galactic disk, 
in this case that of NGC 1058. (see Fig. 1) 

This suggests to me that the velocity dispersion 
of the gas must be determined by physics more or 
less independent of star formation. Star formation 
might then be simply the response of the turbulent 
gas disk to gravity, depending on the amount of gas 
available from smooth or lumpy accretion. 

What is the physical mechanism driving the tur- 
bulence? Mac Low (1999, 2003) estimates the dissi- 
pation rate for a supersonic turbulent fiow to be 
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where p is the average density, wrms the rms velocity, 
and A_D the efl[ective driving scale of the turbulence. 
Note in particular the cubic dependence of the re- 
quired energy input on the rms velocity. A factor of 
two reduction in the velocity will lead to an order of 
magnitude lower energy requirement. 

Sellwood & Balbus (1999) suggested that the 
magnetorotational instability (MRI) (Balbus & 
Hawley 1991, 1998) can drive significant turbulence 



in a differentially rotating galactic disk. They esti- 
mated the energy input by noting that the energy 
input from the MRI 



e = Trd,^, 
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where T^^ is the Maxwell stress tensor, and the 
angular velocity. Numerical models by Hawley et al. 
(1995) suggest that T^^ = O.GB^/Stt, so 

e = (1.2 X 10-2* erg cm^^ s"^) x (6) 
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We can solve equation (4) for the rms velocity that 
could be powered by this level of energy input, get- 
ting wrms ~ 4 km s"^, just a little under the typical 
values observed by Petric & Rupen (2007). 

Several groups have now simulated the MRI in 
galactic disks. Dziourkevitch ct al. (2004) did three- 
dimensional global models with an adiabatic equa- 
tion of state, and found velocity dispersions vary- 
ing with height from 1.5 km in the midplane to 
3 km s"^ at 500 pc and above. Kim et al. (2003); Pio- 
ntek & Ostriker (2004, 2005) and Piontek & Ostriker 
(2007) performed a series of models using a shearing 
box in which they moved from an adiabatic equa- 
tion of state to the introduction of thermal instabil- 
ity to the inclusion of both thermal instability and 
vertical stratification. They concluded that the par- 
tition of the medium into warm and cold phases was 
vital to maintaining significant velocity dispersions, 
with the velocity dispersion in the warm gas ranging 
from 5 km at the midplane to over 10 km at 
500 pc, as shown in Figure 2. 

The MRI seems to provide a floor to the possible 
velocity dispersion in galactic disks. This may ex- 
plain the success of models that follow gravitational 
instability in disks of stars and isothermal gas (em- 
bedded in dark matter halos), which can reproduce 
the Schmidt Law (Li et al. 2005a, 2006), and show a 
tight exponential correlation between the minimum 
value of the gravitational instability parameter Q^g 
(Rafikov 2001) and the star formation timescale (Li 
et al. 2005b). 

3. MAGNETIC REGULATION 

Aside from driving the MRI, magnetic fields can 
also directly affect the location and perhaps the rate 
of star formation in galaxies. Dobbs & Bonnell 
(2007) computed models without self-gravity of mul- 
tiphase gas passing through a spiral potential with 
and without magnetic fields. They found that, with- 
out magnetic fields, clear spurs formed along their 
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Fig. 1. Distribution of velocity dispersion of H I in the galaxy NGC 1058, overlaid on the Ha emission from Ferguson et 
al. (1998) in greyscale. This demonstrates the surprising uniformity and lack of correlation with star formation of the 
velocity dispersion. The regions of highest dispersion are labeled N, C, and S. The x and y axis are the RA and Dec in 
B1950 coordinates. The contours are in km s~^ and start in steps of 0.5 km s~^. Black is used for dispersions between 
5.5 to 7, cyan for 7.5 to 9, green for 9.5 to 11, red for 11.5 to 13, and magenta for 13.5 to 15 km s^^. From Petric & 
Rupen (2007). 
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Fig. 2. Left: Density in a stratified, two-phase, magnetorotationally unstable disk, computed in a shearing- sheet 
approximation. Right: Vertical profile of mass-weighted velocity dispersion Sv for this model for gas within and between 
the stable thermal phases, as well as for the total gas. From Piontek & Ostriker (2007). 



spiral arms, while with fields having magnetic pres- 
sure exceeding the thermal pressure distinct spurs 
were no longer formed, though some interarm struc- 
ture remained. 

Kim & Ostriker (2006) performed shearing box 
models of self-gravitating, magnetized gas interact- 
ing with the fixed stellar potential of a single spiral 
arm, an approximation appropriate to a disk dom- 
inated by stellar mass such as that of the Milky 
Way. They demonstrated spur formation and grav- 
itational collapse as natural consequences of the in- 
teraction of fields with self-gravity, as shown in Fig- 
ure 4. Shetty & Ostriker (2006) followed up on these 
models with global two-dimensional models of self- 
gravitating, magnetized gas in a disk with (and with- 
out) an imposed stellar spiral potential, demonstrat- 
ing the formation of spurs and self-gravitating ob- 
jects that look likely to form giant molecular clouds 
and ultimately stellar clusters. 

In a magnetized disk, collapse can only occur in 
regions where the gas accumulation length along field 
lines is shorter than the generalized Toomre instabil- 
ity length in that region. Cloud formation by such 
collapse will then lead to the surrounding gas having 
reduced mass-to-flux ratio, and thus resisting col- 
lapse more strongly. This feedback mechanism can 
be thought of as magnetic regulation of star forma- 



tion (Shu et al. 2007). 

In spiral galaxies, gas accumulation occurs in spi- 
ral arms, so the work of Kim & Ostriker (2006) and 
Shetty & Ostriker (2006) can be used to derive a 
quantitative estimate of star formation in a magnet- 
ically regulated galaxy. Shu et al. (2007) demon- 
strate that they can derive a Schmidt Law depen- 
dence of the star formation rate on gas surface den- 
sity from considering spiral arms as example accu- 
mulation mechanisms. However, this is a special 
case for the general question of star formation in 
galaxies. It remains to be seen whether this mech- 
anism can explain the behavior of irregular galaxies 
without prominent spiral arms, starburst galaxies, 
or high-redshift, gas-rich objects. 

4. GALAXY FORMATION 

Recently the effect of magnetic fields during the 
galaxy formation era has been directly studied in a 
preliminary adaptive mesh refinement simulation by 
Wang & Abel (2008). They studied the coUapse of 
an isolated halo of total mass 10^° Mq initially hav- 
ing an NFW profile (Navarro, et al. 1996). In this 
first model, further collapse was prevented when the 
Jeans scale approached the grid scale Aa; = 26 pc, 
so no mass was transferred to a collisionless popu- 
lation of stars during the evolution of the galaxy. 
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Fig. 3. Left: Star formation timescale rgp correlates exponentially with the initial disk instability Qsg for both low 
temperature (open) and high temperature (filled) models. The solid line is a least-squares fit. From Li et al. (2005b). 
Right: A comparison of the global Schmidt laws between simulations from Li et al. (2006) and observations. The red 
line is the least-square fit to the total gas surface density of the simulated models, the black solid line is the best fit of 
observations from Kennicutt (1998), while the black dotted lines indicate the observational uncertainty. The color of 
the symbol indicates the rotational velocity for each model; labels from M-1 to M-4 are sub- models with increasing gas 
fraction; and open and filled symbols represent low and high temperature models, respectively. 
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Fig. 4. Surface density (logE/Eo in grey scale) at times of 5.6, 6, and 6.3 orbital times in the magnetized, 3D shearing 
sheet model of Kim & Ostriker (2006). The scale length for typical parameters is ~ Ly/2 — 3.14 kpc. 



As a result, gravitational stirring of the disk is effi- 
cient, maintaining velocity dispersions of 10 km s~^ 
or more, but not from an astrophysically relevant 
mechanism. Magnetic field was assumed to initially 
thread the halo with strength of a nanogauss, under 
the assumption of field production by previous gen- 
erations of star formation and active galactic nuclei 
(e.g. Rees 2006). 

The initial evolution of the galaxy shows rapidly 
increasing amounts of gravitationally unstable gas 
that ought to produce large numbers of stars in the 
first several hundred Myr. During this period, the 
field is amplified in the disk, reaching a few per- 
cent of equipartition after about 500 Myr. The field 
only begins to play a role in supporting the gas af- 
ter that point, however, while star formation pro- 
ceeds with little difference between pure hydrody- 
namic and magnetized models for the entire previ- 
ous period. The initial conclusion appears to be that 
even with relatively large initial fields, magnetic ef- 
fects can be neglected during the initial period of 
galaxy formation, until a large-scale field has been 
generated. 

5. FIELD GENERATION 

This brings us to the question of how such 
large-scale magnetic fields are formed. Classical 
Alpha-Omega dynamos require many Gyr to gener- 



ate observed fields, although somewhat shorter times 
may be given by a cosmic-ray driven dynamo (see 
Otmianowska-Mazur, this volume, electronic edi- 
tion). Turbulent dynamos can generate strong fields 
much more quickly, but the dominant scale for field 
structure is then the diffusive scale (Schekochihin et 
al. 2004), so there is no large-scale, coherent field 
structure like that observed in galaxies. However, 
the question of how the tangled field generated by a 
turbulent dynamo evolves when embedded in a dif- 
ferentially rotating galaxy remains unaddressed. 

The simulations I have discussed in this paper do 
not answer this question, but they do offer sugges- 
tive results that might point the way towards the an- 
swer. Beginning with the simplest case, we note that 
the in their stratified, shearing box MRI simulation 
with thermal instability of Piontck & Ostriker (2007) 
shows field growth to 3-4 /iG in the cold gas within 
300 Myr, and then maintains that field strength for 
another Gyr. Nishikori et al. (2006) performed a 
three-dimensional, resistive MHD model of a global 
torus with initially weak azimuthal field threading 
10^ K gas. They found an increase in magnetic en- 
ergy of eight orders of magnitude over a period of 
200 Myr, and then maintain the field at that level 
for another 600 Gyr. 

Wang & Abel (2008) find a similar behavior for 
the magnetic energy in their simulation of disk for- 
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Fig. 5. Left: Time history from the AMR model of Wang & Abel (2008) of disk gas kinetic energy (squares), thermal 
energy (diamonds), and magnetic energy (asterisks) , showing the strong growth of magnetic energy over time. Right: 
Horizontal and vertical slices of magnetic pressure after 1.088 Gyr of evolution. The plot square has a side of 11 kpc. 
From Wang & Abel (2008). 



mation in a collapsing halo, showing a factor of 10^ 
amplification of field energy over 600 Myr, as shown 
in Figure 5. They further demonstrate that the mag- 
netic field in their disk indeed has large-scale struc- 
ture after 600 Myr, with several field reversals, clear 
spiral structure, and field strengths in spiral arms 
exceeding 10-20 ^J.G (Fig. 5). 

6. SUMMARY 

The major points I have made in this paper 
arc that the MRI may offer a mechanism to main- 
tain a minimum level of turbulence regardless of the 
strength of star formation. This might then sup- 
port the picture that gravitational instability acting 
on disks regulates the strength of star formation in 
galaxies. The alternative view that magnetic fields 
regulate star formation appears likely to be valid at 
least morphologically, but it remains to be seen if the 
rate of star formation requires magnetic regulation 
to explain. 

A related question is how the large-scale mag- 
netic field in galaxies is generated. Small-scale tur- 
bulent dynamos to generate field strength, perhaps 
with large-scale coherent fields produced at the end 
by classical dynamo action is an intriguing specu- 
lation that remains consistent with the increasingly 
detailed models that I have reviewed here. 
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